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ON GLOBAL NON-OSCILLATION OF LINEAR ORDINARY 
DIFFERENTIAL EQUATIONS WITH POLYNOMIAL 
COEFFICIENTS 

DMITRY NOVIKOV AND BORIS SHAPIRO 


Abstract. In this note we show that a linear ordinary differential equation 
with polynomial coefficients is globally non-oscillating in if and only if 

it is Fuchsian, and at every its singular point any two distinct characteristic 
exponents have distinct real parts. As a byproduct of our study, we obtain a 
new explicit upper bound for the number of zeros of exponential polynomials 
in a horizontal strip. 


1. Introduction 

Let us recall the classical notions of disconjugacy and non-oscillation of a linear 
ordinary differential equation, see e.g. [3]. 

Definition 1. A linear ordinary differential equation of order k 

ak{z)y^'^^ + + ... + aQ{z)y = 0 , ( 1 ) 

with continuous coefficients aj{z), j = 0,... ,k defined in a neighborhood of some 
simply-connected subset I of R or C, is called disconjugate (resp. non-oscillating) 
in J, if every its nontrivial solution has in I at most k—1 zeros (resp. finitely many 
zeros) counted with multiplicities. 

Observe that every equation o is disconjugate in any sufficiently small interval 
in R (resp. any sufficiently small disk in C) centered at an arbitrary point zq G ®. 
(resp. zq G C) such that ak(zo) ^ 0. Analogously, every equation ([T]) is non¬ 
oscillating in any compact simply-connected set free from the roots of ak{z). 

The study of different aspects and criteria of disconjugacy and non-oscillation 
has been an active topic in the past. While there exist satisfactory criteria of 
disconjugacy for the second order equations, the situation with the higher order 
equations is more complicated. A number of necessary/sufficient conditions of 
disconjugacy for subsets of M and C are known in the literature mostly dating 
back at least four decades, see e.g. m, 0,0- In the case of equations of order 
2, disconjugacy is closely related to Sturm separation theorems; for higher order 
equations there is a related version of multiplicative Sturmian theory developed in 

[la¬ 
in this paper, for a linear differential equation with polynomial coefficients, we 
introduce the notion of its global non-oscillation in CP^ by which we mean its 
classical non-oscillation in an arbitrary open contractible domain obtained after 
the removal from CP^ of an appropriate cut connecting all the singular points. 
Although oscillation/non-oscillation in the complex domain have been studied since 
the 1920’s, (see e.g. 0), the notion of global non-oscillation seems to be new. As 
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an experienced reader can easily guess, the main motivation for our consideration 
comes from the second part of Hilbert’s 16th problem. 

Consider a linear homogeneous differential equation 

+ ... + Poiz)y = 0 , ( 2 ) 

with polynomial coefficients Pk{z), Pk-i{z ),..., Po{z), and GCD{Pk, Pk-i, ■ ■ ■, Po) = 
1. Let S be the set of all singular points of (HI) in CP^, i.e., the set of all roots of 
Pk{z) (together with oo if some of the limits lim^^oo z^ Pk-j{z)/Pk{z), j — 0,k 
is infinite). For a given equation (HJ, let d denote the cardinality of S. 

Definition 2. A system C := {Cj}jZl of smooth Jordan curves in CP^, each of 
them connecting a pair of distinct singular points, is called an admissible cut for 
equation HD if and only if: a) for any i ^ j, the intersection Ci C\Cj is either empty 
or consists of their common endpoint; b) the union is topologically a tree in 

CP^, i.e., the complement CP^ \ UjC^ is contractible; c) each Cj has a well-defined 
tangent vector at each of its two endpoints. 

In particular, there exist admissible cuts consisting of straight segments connect¬ 
ing the singular points of ©• 

Definition 3. Equation (HD is called globally non-oscillating if, for any its admis¬ 
sible cut C, every its nontrivial solution has finitely many zeros in C\C. 

The main result of this paper is the following criterion of global non-oscillation. 

Theorem 4. Equation HD is globally non-oscillating if and only if: 

(i) it is Fuchsian; 

(ii) at each singular point all distinct characteristic exponents have pairwise 
distinct real parts. 

Remark 5. One can easily notice that HD is globally non-oscillating if and only 
if some (and therefore any) domain CP^ \ C can be covered by finitely many open 
disconjugacy domains. Observe that if one knows such a covering, then one gets an 
immediate upper bound for the total number of zeros of nontrivial solutions of (HD 
in CP^ \ C. Namely, if the number of open disconjugacy domains covering CP^ \ C 
equals I, then any nontrivial solution of (HD has there at most {k—l)l zeros counted 
with multiplicitties. 

In view of Remark HI the following problem is of fundamental importance. 

Main Problem. Given an arbitrary equation HD satisfying the assumptions of 
Theorem HI estimate from above the number of disconjugacy domains which can 
form an open covering of CP^ \ C, for some admissible cut C. 

Observe that in case of a Schrbdinger equation 

-y” + P{z)y = 0 

with a polynomial potential P{z), there is a classical construction of such coverings 
using the Schwarzian derivative of two linearly independent solutions of the latter 
equation which goes back to R. Nevanlinna, m- 

Remark 6. Let us also mention that Proposition HI below, which is an important 
technical tool used to prove Theorem HI is a new result in the classical area of the 
upper bounds for the number of zeros of exponential polynomials and, therefore, 
it is of independent interest. Such upper bounds are required in a wide range of 
mathematical disciplines, from applied mathematics to number theory. Essential 
progress in this area has been made in the 70’s in the papers Ha, [E], [S]. But, to 
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the best of our knowledge, in all the previous literature one only considered compact 
subdomains in C, mainly disks and rectangles, while Proposition IH] considers the 
case of an infinite strip. 

Acknowledgements. The second author is grateful to the Department of Mathemat¬ 
ics and Computer Science of the Weizmann Institute of Science for the hospitality 
in January 2010 and February 2015 when this project was initiated and carried out. 
The first author wants to thank G. Binyamini for many discussions of the upper 
bounds of the number of zeros of Fuchsian equations over the years. 

2. Proofs 

Our proof of Theorem 0] consists of several steps. 

Step 1. The necessity of Conditions (i) and (ii) for global non-oscillation of equa¬ 
tion ( 01 ). 

Indeed, if m has a non-Fuchsian singularity at p G CP^, then, for any suffi¬ 
ciently small e > 0 , almost any solution of (0]) has infinitely many zeros in the 
e-neighborhood of p with a removed straight segment connecting p with some point 
on the bounding circle. This property contradicts to global non-oscillation. To 
finish Step 1, consider a Fuchsian singularity of (0]) with two distinct characteristic 
exponents of the form a -I- bil and a + b 2 l. Then there exists a solution of (0]) with 
the leading term 2 cos { ^2 ■ Such a solution has infinitely many ze¬ 

ros accumulating to p which are located close to the horizontal line passing through 
p. This again contradicts to global non-oscillation. 

Step 2. Reduction to small neighborhoods of singular points. 

For any sufficiently small e > 0, construct a simply-connected domain Ue C CP^ 
by: a) taking the large disk {jzl < e~^} with the e-neighborhoods of all zeros of 
Pk removed, b) making cuts by straight segments between the bounding circles so 
that the obtained domain becomes contractible. 

The following complex analogue of the classical de la Vallee Poussin theorem [1] 
is proved in [T71 Theorem 2.6, Corollary 2.7]. 

Lemma 7. Consider a homogeneous monic linear ordinary differential eguation 
with holomorphic coefficients 

+ Ofe_i(t) H-h ao{t) y = 0, t G C. 

Then the variation of the argument of any solution y{f) along a circular arc 7 of 
a known length is explicitly bounded in terms of the uniform upper bounds Ai = 
sup(g^ |a*(t)|, i = 0,...,k-l. 

Lemma [7] implies an explicit upper bound B(e) for the number of zeros of any 
solution of (01 in (7^. More exactly, the upper bound will depend on the upper 
bounds on the restrictions of aj to dU^. The latter are polynomial in e~^ if the 
coefficients aj are polynomials, so the upper bound is also polynomial in e~^. 

Remark 8. Observe that, for any admissible system of cuts C and any sufficiently 
small e, the domain CP^ \ C can be covered by finitely many Ue (choosing different 
straight lines connecting the bounding circles) and finitely many sectors of finite 
radii centered at the singular points of (0]). This observation reduces the proof of 
Theorem 0] to providing finite upper bounds for the number of zeros of solutions of 
@ in these sectors, see below. 

Step 3. Equations with constant coefficients. (’’Reduction" to the case of equations 
with constant coefficients in a neighborhood of a Fuchsian singularity is obtained 
by using the logarithmic chart centered at the singularity. See also Steps 4-5.) 
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Proposition 9. For any a > 0 and for any equation 

EQ : + ... + a^y = 0, Oj G C, afe 7 ^ 0 (3) 

such that all its distinct characteristic roots have distinct real parts, 

(1) there exists an upper bound i{EQ, a) for the number of zeros of all nontrivial 
solutions of in the horizontal strip {IIq : |3(z) < a|}. (Here zeros are 
counted with multiplicities.) 

(2) in the generic case when all roots Xj, j = of the characteristic 

equation of ( 0 ^ are simple, with JRAi < KAi < • • • < KAfc, we get 

UEQ, a)<{k- If + -(fc - 1)£(EQ) [a(S + 2) + 0 log 4], (4) 

TT 

where £(EQ) is the length of the shortest polygonal path passing through all 
Xj and 

0 ■= “ — , max 

The case of multiple characteristic roots will be considered in Step 4. 

Our approach to the proof of Proposition |9] is inspired by the Wiman-Valiron 
theory, see M- The main construction below has a strong resemblance with the 
notion of a tropical polynomial in the modern tropical geometry. The proof itself 
is based on a rather long sequence of lemmas and the last argument is given at the 
end of Step 5. 

The general solution of is given by: 

y = ^ Aj{z)e^^^, where deg Aj{z) = nj,''^{nj + 1) = fc. (5) 

3 

Define the domain of a single term y-dominance in fla as 
G{y,a) :={zGn„ |3j=j(^),3e>0 : > (1-e) ^ |Al,(z)e"^^|}. ( 6 ) 

ii^3 

Note that G{y,a) may contain at most minnj < k zeros of y, namely the 
common zeros of all Aj{z). In particular, G{y,a) contains no zeros of y at all in 
the case of simple characteristic exponents. 

Lemma 10. The complement IIq, \ G{y,a) can be covered by at most k -\- kf' + k^ 
horizontal boxes (of height 2a) of the total width not exceeding 

k^{k + 1)(40 InA: + daS + 4a) + 8k^Q. 

We first consider the case of simple characteristic exponents Xj. This case is more 
transparent and the resulting estimates seem to be of correct order of magnitude. 
In this case the polynomials Aj{z) are constants and will be denoted by aj. 

Lemma 11. In the case of simple characteristic exponents Xj, the complement 
IIq \ G{y, a) can be covered by at most k — \ horizontal boxes (of height 2a) of the 
total width not exceeding 

2a(fc-I)S + 2(/c-l)01n4. (7) 

The principal case in Lemma [TT] is a = 0, i.e. Hg = ffi.. 

Lemma 12. In the above notations, R \ G{y, 0) is contained in the union of at 
most k — 1 closed intervals of the total length less than or equal to 2{k — 1 ) In4 • 0 . 


’AXj — SlAj+i 
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To prove Lemma m we need an additional statement. In with coordinates 
consider the 1-parameter family {Ptj{u)}^_^ of k points given by 

fj, = K(Aj), (j) = \w |aje^^”|, 

where it is a real-valued parameter. For a given value of it G K, introduce the 
piecewise-linear function as the least concave majorant of By 

this we mean the minimal concave function (jjuilJ-) defined in the interval [5i(Ai), 3i(Afc)] 
such that all points {Ptj{u)}’j^^ lie non-strictly below its graph, i.e. have their (p- 
coordinate smaller than or equal to that of (One can easily see that the 

graph of (puil^) is the upper part of the boundary of the convex hull of {Ptj{u)}j^i 
connecting Pti{u) and Ptk{u).) Observe that, for any it G M, 

^„(/i) = 0 o(m) + MM- (8) 

Lemma 13. If, for j = 1,..., fc — 1, 

|(^„(5R(A,+i))-<(.„(»(A,))| >ln4, 

then u G G[y, 0). 

Proof of Lemma \PA Define the central index of (f>u (m) by the formula: 

i(it) := {i |5ft(Ai) is the point of the global maximum for 
comp. Ch. 1, [14] . Then for any j ^ i, 

< exp(</.„(KA,)) < 4-b-*(“)l 

Therefore the inequality in the definition ([5]) of G{y, 0) follows after the summation 
of a geometric series. □ 


Corollary 14. If —u lies outside the In 4 • Q-neighborhood of the set of all slopes 
of <Po{pl), then u G G{y,0). 

Proof. Formula (I5|) implies that each slope of pu (m) equals the sum of the respective 
slope of <po{y) and it. Therefore in the considered case, the absolute values of 
all slopes of (puih-) exceed In4 ■ 0 , and the statement follows immediately from 
Lemma M □ 


Proof of Lemma\TB. The In 4 • 0-neighborhood of the set of slopes of (poiff) consists 
of the union of at most fc—1 intervals of total length not exceeding 2(fc—1) In 4-0. □ 


Proof of Lemma \Tl\ Consider the general case of Lemma [TT] with a > 0. 

We repeat the above construction of Lemma[T2]for z running along the horizontal 
line 3z = V with |i;| < a. For every fixed v, consider in with coordinates (/i, p), 
the 1 -parameter family {Pt'j{u)}j^i of k points given by 

/i = {K(A,),(^ = ln|aje^^(“+^'')|} 


where u is a real parameter. Introduce Puih) the least concave majorant of 
{Pt'-{u)}^^i, for a given value of u G R. Observe that, for any it G R, 

+mm- (9) 

Now consider the set Sla '■= ^-a<v<a{kj{v)}, where kj{v) are the slopes of 
Poifi). We claim that Sla is the union of at most k — 1 closed intervals. Indeed, 
the set of slopes {kj{v)} changes continuously with v, and consists of no more than 
k — 1 points for each fixed v. 

Moreover, as In = In ja^l — v'HXj, the points {Pt'j{0)}j^i defining Pq{pl) 

depend linearly on v, namely they move up or down as v changes. The inequality 


dkj{v) 

dv 


< E 
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is straightforward. Therefore the total length of Sla is at most 2a'B.{k — 1). 

By Corollary [lU if —u lies outside the In 4 • ©-neighborhood of Sla, then, for 
any |r;| < a, u + Iv lies in G{y, a) which settles Lemma fTTl □ 


Step 4. Case of multiple characteristic exponents. 

In this case the dependence on v of (analogs of) points Pt'-{u) seems to be more 
complicated, and we are forced to consider the slopes of all chords connecting these 
points, and not only those which lie on the boundary of their convex hull. This 
apparently leads to an excessive upper bound of the total width of IIq, \ G{y, a). 


Proof of Lemma \1(A Consider the absolute value of the ratio of any two terms 
in ®. The complement Ha \ G{y, a) lies in the union S of the sets 11°^/ = 
{\\nrjji{z)\ < Infc}, where rjji is the absolute value of the ratio of two terms 
in (O. 


We can write 

where 

Set W = {|lR(z 
we have 

Additionally, 




( 10 ) 


^33' ~ ^jj' ~ 

Zi)\ > 4fc0}, where Zi runs over all roots of all Aj. Outside W 


A 

du 


\n\Aj/Aj, 


dv 


In \Aj/Aj, 


< 


\ 03 . 


= {u + Iv e n„ : |lnrjy(u)| < Ink + + a\9jj,\}, 

outside W. Note that Sjy is the union of boxes, since its definition is independent 
of V. 

Therefore, 

du 

outside W. 

Thus 'Sjji intersects each connected component of R \ W in an interval of length 
at most 4:\9jji In fc -b 4a|^jj/1 + 4a. In other words, \ IT is the union of at 
most fc -b 1 boxes of total width not exceeding (fc -b 1)(4 In k\9jjr |-b 4a|^jy' | -b 4a). 

Taking the union over all possible pairs (j, j'), we conclude that S \ IT lies in the 
union of at most fc^(/c-bl) boxes of total width at most fc^(fc-bl)(40 In fc-b4aS-b4a). 
As IT n IIq is the union of at most k boxes of width at most 8fc0 each, we obtain 
that S lies in the union of at most fc -b -b boxes of total width at most 
fc^(fc-b 1)(40 Infc-b 4aS-b 4a)-b 8fc^0. □ 


> 


\^33 


( 11 ) 


Finally let us explain how Lemmas (TU] and [TT] imply Proposition [5] Consider the 
space QPa = ^ *11^} oldimension fc = X)(l + 'l6g^j) consisting 

of exponential polynomials, where A = {Aj} C C is some finite set. The following 
result was proven in [B]. 

Theorem 15 ([B]). The number of zeros of any function f € QPa in n bounded 
convex domain U does not exceed 

fc-I + -/:(A)diam(t7), (12) 

TT 

where T(A) is the length of a shortest polygonal path passing through all points of 

A. 
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Theorem [TS] immediately implies an estimate on the number of zeros of y in 
the boxes Bj of Lemma uni and HH In the case of simple characteristic exponents 
(second part of Proposition 19]) 

diam.Bj < 2{k — 1) [aS + 0 log 4] + 4(A: — 1)q;, 

and O follows. 

Step 5. Equation with non-constant coefficients in a semistrip. 

In general, solutions of ([2]) considered in the logarithmic chart near its Fuchsian 
singularity have the form 

2/= (13) 


— X/ ’’(1 + ej,r), 

r—0 

and Cj^r is 27r/-periodic, ej^r = O(e^) in any semistrip = {|3z| < a, < /3} 
for some f3 depending on m only. To simplify our notation, let us assume that 
/3 <0. 

Lemma 16. Assume that \ej^r\ < in Let Aj{z) = Y^Lo0.3,rz'^^ ’’ and 
W he as in the proof of Lemma [7^ Then 

I log \Aj/Aj\ < CEQ,a in \ W, (14) 

where CEQ,a is some constant depending on a and ([7]) only. 

Proof. Let A{z) = \aj,r\\z\'^- Evidently, Aj{z) < + \zm\), where Zm are 

the roots of Aj{z). Also, \Aj — Aj\ < Ce^’^Aj{z). Therefore, for z G \ W, we 
get 


\log\Aj/Aj \ < Ce 


Vtz 


A(z) 

\M^)\ 


< Ce^^ ■ 


z - z„ 


< 


(15) 




n 

|Zm|<2|2| 


I + \z„ 


z - Z„ 


n 


|z/z„ 


II 


\zm,\>2\z\ 


1 - |z/z„ 




\4ke) 


k-t 


Clearly, the latter function can be majorized by some number depending on (7, a 
and k, 0 only. The constants C,k,Q are determined by ©. □ 


Remark 17. Actually, dependence of CEQ,a on a is very simple (as 0{a’^) as 
a ^ oo), but we do not need this. 

Lemma 18. In the above notation, the zeros ofy in IIo.^^ lie in at most k + k^ + k^ 
boxes of total width at most 

k^{k + 1)(40 In fc + laS + 4q; + iCEQ,a) + 8fc^0. 

Proof. We repeat the proof of Lemma (TUI Namely, consider the absolute value 
fjji of the ratio of any two terms in (1131) . The complement IIq; \ G{y,a) lies in 
the union E of the sets T,°j, = {| In f n'iz)\ < Infc}. But, according to Lemma |16l 
I logrjji — logrjji I < CQE,a, where rjji was defined in the proof of Lemma fTHl So, 
it is enough to require | \a.fjji{z)\ < In/c -|- CEQ,a, i-e. outside W 

C Ejj' = {m + /u G IIq : |lnrjy/(M)| < Ink + a\f,jj'0jj' \ +a\9jy\ +CEQ,a}- 

Repeating the same arguments as in Lemma [TU] with 'Ejj/ instead of 'Pjj', we arrive 
at the required estimates. □ 
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Proof of Theorem\^ Let + ... + bky = 0 be the reduced form 

(=divided by its leading term) of (HJ in the logarithmic chart near its Fuchsian 
singularity. Assume that bj{z) are bounded by C in (The Fuchsian property 
implies that each bj{z) tends to some finite limit when z ^ oo in 

Example in HZ] immediately following after Corollary 2.7 of this paper, implies 
that y has at most 2 {k + 1) + zeros in where 

(. < 2k^{k + 1)(40Infc + 4a5 + 4a + 4CEQ,a) + 16fc^0 + 4{k + + k^)a, 

is the total perimeter of all boxes appearing in Lemma [T 8 | 

After going back from logarithmic chart to the original coordinate, we obtain an 
upper bound for the number of zeros of any solution of ([5]) in the sector {|z — p| < 
e^, I argz| < a} at the Fuchsian singular point p. 

□ 


The sequence of steps 1-5 settles Theorem |4| 
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